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We explicitate the relation between Hamiltonians for networks of interacting qubits in the XY Z 
model and graph Laplacians. We then study evolution in networks in which all sites can communicate 
with each other. These are modeled by the complete graph K n and called all-to-all networks. It 
turns out that K n does not exhibit perfect state transfer (PST). However, we prove that deleting 
an edge in K n allows PST between the two non-adjacent sites, when n is a multiple of four. An 
application is routing a qubit over n different sites, by switching off the link between the sites that we 
wish to put in communication. Additionally, we observe that, in certain cases, the unitary inducing 
evolution in K n is equivalent to the Grover operator. 



I. INTRODUCTION 

Networks of interacting qubits are a generalization of 
spin chains. These are of theoretical importance for the 
study of many-body quantum systems and could consti- 
tute a good test ground for technologies spanning from 
quantum key distribution in multi-user networks to var- 
ious nano-scale devices. 

Since the first works @, [H, (see also @ , for a re- 
view), networks of interacting qubits are considered to be 
good candidates for engineering perfect quantum chan- 
nels and allow information transfer between distant par- 
ticles. In this perspective, such networks appear to be 
useful for the implementation of data buses in quantum 
mechanical devices, in particular because undergoing a 
free dynamics after an initial set-up. One of the prob- 
lems arising in the scenario is given by natural disper- 
sion effects, which determine a loss of information often 
proportional to the distance (i.e., the number of spins) 
between communicating sites. Ways to circumvent the 
issue are based on a local tuning of the couplings [5| or 
protocols for lifting the encoding of qubits into multipar- 
ticle states 0, 0| (see also @ and the references therein) . 

Another approach is to study what network topologies 
guarantee a small decay of fidelity, when the couplings 
are homogeneous and constant during the entire evolu- 
tion of the system. Results in this direction have pro- 
vided examples of the so-called perfect state transfer (for 
short, PST), in relation to combinatorial properties of 
the graphs modeling the networks (a list of references on 
this area is in [l[). In the XY model, when considering 
a single excitation, it has been shown that PST depends 
essentially on the eigensystem of the adjacency matrix of 
the graph, because certain invariant eigenspaces of the 



total Hilbert space evolve independently. In the cases 
analyzed so far, even if PST occurs between two specific 
sites of a network, routing arbitrarily over the entire net- 
work still requires an external controller. 

Here we will deal with PST in certain networks of 
spin-half particles with the XYZ interaction. Specifi- 
cally, we will describe the structure of the Hamiltonian 
which governs the single excitation setting. This turns 
out to be proportional to the (combinatorial) Laplacian 
of the graph modeling the network. Thus, it may be 
worth to remark that Hamiltonian arising in the XY and 
theX YZ interactions are naturally associated to the ad- 
jacency matrix and the Laplacian respectively and that 
these objects are the two most common matrix represen- 
tations of graphs. 

We will study quantum evolution in a network in which 
all sites can communicate between each other (in both 
directions). Such a network is modeled by the complete 
graph on n vertices K n and it is called all-to-all network. 
The number of links in K n is n(n — l)/2. It turns out 
that K n does not exhibit PST. However, we will show 
that deleting an edge in K n will allow PST between the 
two non- adjacent sites in the obtained network, when n 
is a multiple of 4 (apart the trivial case n = 2). This is 
modeled by a graph usually denoted by K~ , an all-to- 
all network with a missing link. The phenomenon is due 
to interference effects and it is counterintuitive. On the 
side, we remark that in certain cases the unitary inducing 
evolution in K n is the Grover matrix up to an overall 
phase. 

An application is routing a qubit on K~ (i.e., trans- 
ferring the qubit between any two vertices of K~ , arbi- 
trarily). The protocol involves an external agent, whose 
role is to switch OFF only the link between the sites that 
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we wish to put in communication. When the process in- 
volves all n sites, it requires exactly 2n — 1 ON/OFF 
switching operations. This does not provide any direct 
advantage over considering the empty graph as a network 
(this is trivially the graph with no edges), with a single 
edge ij if communication between site i and site j is de- 
sired. On the other hand, dynamics on K~ is far more 
complex and the same edge configuration is potentially 
useful for other tasks beyond perfect communication of 
a qubit. 

The structure of the paper is as follows. In Section 2 we 
give minimal background information on the XY Z model 
and the necessary mathematical definitions. We will clar- 
ify the exact relation between the XY Z Hamiltonian and 
the Laplacian matrix (this relation was mentioned in Q , 
but not given explicitly) . In Section 3 we study evolution 
on K n . Section 4 is devoted to K~ . Wi will also mention 
a generalization in which we delete an arbitrary number 
of vertex disjoint edges. A brief conclusion is drawn in 
Section 5. 



II. SET-UP 

Let G = (V, E) be a simple undirected graph (that is, 
without loops or parallel edges), with set of vertices V(G) 
and set of edges E(G). We take V(G) — {1, ...,n} and 
assume that |-E(G)| = m. The degree d(i) of a vertex i is 
the number of edges incident with i. The adjacency ma- 
trix of G is denoted by A(G) and defined by [A{G)]ij = 1, 
if ij £ E(G); 0, [A(G)]ij = if ij $ E(G). 

The adjacency matrix is a useful tool to describe a 
network of n spin-half quantum particles. The parti- 
cles are usually attached to the vertices of G, while the 
edges of G represent their allowed couplings. If one 
considers the XYZ interaction model (isotropic Heisen- 
berg model), then {i,j} £ E(G) means that the par- 
ticles i and j interact by the Hamiltonian [H(G)]ij — 
{XiXj + YiYj + ZiZj), where Xi, Yi and Zi are the Pauli 
operators of the i-th particle (here we consider unit cou- 
pling constant). In extension of the XY model, we decide 
to call this the XYZ model, as the Z interaction has been 
added to the XY Hamiltonian. Thus, the Hamiltonian 
of the whole network reads 

Hxyz(G) = X - Y^[A(G)]i3 (XiXj + YiYj + Z^) , (1) 

and it acts on the Hilbert space /C = (C 2 )®". 

Let us now restrict our attention to the single exci- 
tation subspace H = C™, i.e., the subspace of dimen- 
sion n spanned by the vectors {|1), . . . , A vector 
|j) indicates the presence of the excitation on the j-th 
site and the absence on all the others. This is equiva- 
lent to the following tensor product of the Z eigenstates 
I P . . . 010 . . . ) , being 1 in the j-th position. In the basis 

n 

{|1), . . . , \n)}, the Hamiltonian coming from Eq. ([1} has 



the following entries 

[HxY(G)]ij = 2[A(G)]ij, i + j (2) 

and 

[H z (G)]u = i X)[A(G)]y - 2 (3) 

id j 

Eq. © comes from the term X t Xj + Y t Yj and Eq. © 
comes from the term Z^Zj. It is then clear that in the XY 
model the relation between the Hamiltonian and adja- 
cency matrix simply reduces to [HxY{G)]ij = 2[A(G)]ij. 
Here, H X yz(G) = H XY (G) + H Z (G). By associating the 
vertex i € V(G) to the vector \i) 6 {|1), |n)}, we can 
introduce the following modified version of the adjacency 
matrix: 

Definition 1 The XYZ adjacency matrix of a graph G 
is denoted by H(G) and defined by 

f 2, ifijeE(G); 
[H XY z{G)}ij = {^ ifij£E{G)- 
[m- 2d(i), ifi = j. 

We shall drop the subscript XYZ because in this paper 
we only deal the XYZ model. The matrix H(G) has a 
neat relation to the Laplacian of G. Let A(G) be an 
n x n diagonal matrix such that [A(G)]^ = d(i). The 
(combinatorial) Laplacian of G is the matrix L(G) := 
A(G) — A(G). Let I n be the n x n identity matrix. 

Lemma 2 The XYZ adjacency matrix of a graph G is 
H{G) = ml n — 2L(G). 

Proof. By Definition [Q H(G) = 2A(G) + ml n - 2A(G). 
By combining this fact with the definition of Laplacian, 
we obtain the desired expression. ■ 

Corollary 3 Given any graph G, the matrices L(G) and 
H(G) have a common set of eigenvectors. Moreover, if \x 
is an eigenvalue of L(G) then X — m—2fi is an eigenvalue 
of H(G). 

Let {|1), \n)} be the standard basis of an Hilbert 
space H = C\ We associate the vertex i g V(G) to 
the vector \i) £ {|1), \n)}. Let i — V^I. Given a 
graph G, two of its vertices i and j, and a real number 
< t < oo, the fidelity at time t between vertex i and 
vertex j is defined by f G (i,j,t) ■= \ {j\e~ Lli(G)t \i) \. We 
say that two vertices i and j of a graph G admits perfect 
state transfer {w.r.t. the XYZ adjacency matrix) if there 
exists t for which f G {i,j,t) = 1. Let U t {G) = e ~ LHt - G ^ 
be the unitary matrix associated to H{G) as a function 
oft. Therefore, we have f G (i,j,t) = 1, if \[U t (G)]i,j\ = 1. 

III. ALL-TO-ALL NETWORKS 

In this section, we consider a network in which all sites 
interact directly. This network is modeled by the com- 
plete graph K n on n vertices. In K n , we have d(i) = n—1, 
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for every i G V(G). The principal result of the section is 
the following statement: 

Theorem 4 Let K n be the complete graph on n vertices. 
For every vertex i G V(K n ) and value n G N, we have 
the following: 

• mm t f Kn (i,i,t) = 1 - 2/n and t = ^ + for 
k > 0; 

• max ( fjc n (i, i, t) = 1 and t = ■nk/n, for k > 1. 

For every two distinct vertices i,j G y(if„) and waZue 
n G N, we ftove i/ie following: 

• min t /jf n (i, j, t) = is attained when fx n (i.i,t) is 
maximum; 

• maxj /fir n ("*,.?,£) = 2/n is attained when fx n (i 1 i,t) 
is minimum. 

Proof. Since d(i) = n — 1, for every i G V(K n ), we have 
A(iT n ) = J n — I ni where J n is the n x n all-ones matrix. 
Then, L(K n ) = (n - l)J n - ( J„ - J n ) = n/„ - J„. Given 
|S(if n )| = (™), it follows that 

H(Xn) = (fj I n - 2 (nJ„ - J n ) = n ( n ~ 5 h n + 2J n . 

The eigenvalues of H(K n ) are then A^ 1 ' = n(n— l)/2 and 
A^" 1] ([i?(^n)]M - = n (n - 5) /2. Let us de- 

note by J„ the all-ones vector of dimension n. Since 
[H(K n ), J n ] = 0, these two matrices share a common set 
of eigenvectors. The matrix J n is diagonalized by F n , the 
Fourier transform over the group Z„: [F n ]ij = e 1,2 ™ 1 ^ 71 = 
lj 70 . As a consequence, the single eigenvector correspond- 
ing to the eigenvalue Ai is of the form |Ai) = n~ x / 2 J„; 
the n — 1 eigenvectors corresponding to A2 are of the form 
|A|) := n -1 / 2 5Z? = x w^Ij), for every i = 1, n— 1. Now, 
we can read Ut(K n ) = e - iH ( K n)t j n jt s spectral decom- 
position: 



is attained for each t such that |[f/t(-Kra)]ii| is maximum 
and viz. ■ 



The results is that there is no PST in K n (but for the 
trivial cases n ^ 1,2). This is equivalent to say that 
a qubit can not be accurately transferred between two 
sites in a network in which all sites are connected to 
each other. The fidelity given by the evolution in K ni 
with n = 4, 8 is illustrated in the figures below (left and 
right, resp.) for t G [0, ir]. The dashed lines represent 
\[U t (K n )]ij\ with i ^ j; the solid ones, \[U t (K n )] u \. 




For t = zll -f JL j Theorem |4] prompts to the following 
facts. If n is odd then 



[U t (K n )}i 



Sij — 2/,/n, if n = Ak — 1; 
% - 2/n, if n = Ak + 1. 



If n is even, we need to distinguish two cases: 
• n = 4fc, 



[Ct(#»)]« = 



— (Sij — 2/n) y/i, if fc even; 
(oVj — 2/n) y% if fc odd; 




(Sij — 2/n) v% if even; 



(5y — 2/n) y/i, if fe odd. 



U t (K n ) = n- 1 \^e- cXlt J n + e 
Then 

[Ui(^n)]« = n-^- 4 ^)* + 
for every i and 

[U t (K n )] t] = n- 1 (e-<*> 



n — 1 ; n 

-l.X 2 t /.,*J7^* fc 



^ZZJ Ife |j)(fc| 



n - 1 . n(n-5) , 



■"("-5). 



if i 7^ j. Thus, we have 1 - 2/n < \[U t (K n )} u \ < 1, for 
every z. The minimum is attained when t = ^ + 
The maximum is attained for each t = nk/n. So, < 
\[Ut(K n )]ij I < 2/n. For these entries, the minimum value 



This is the Grover operator with an overall phase. Re- 
call that this operator plays a central role in the Grover 
algorithm for database search. The unitary Ut(K n ) in- 
duces a continuous-time quantum walk on K n . Notably 
an alternative version of this search technique has been 
designed as a discrete-time quantum walk algorithm 0. 



IV. ALL-TO-ALL NETWORKS WITH A 
MISSING LINK 

In this section, we consider a network in which all sites 
but two interact directly. This network is modeled by 
K~, the complete graph minus an edge on n vertices. The 
graph K~ is obtained from K n by deleting an arbitrary 
edge. 
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Theorem 5 Let K~ be the complete graph minus an 
edge on n vertices. For every vertex i £ V(K~) and 
value n = 4k, with k £ N, we have the following state- 
ments: 

If i = l,n then 

• min t f K - (i, = and t = f + ^ , for k > 0; 

• max( f K - (i, i, t) = 1 and t = Ttk/2, for k > 1. 
If i ^ 1 , n then 

• mint f K - (i, i, t) = 1 - 2/n and t = ^ + ^, for 
k>0; 

• max ( /g— (i,i,t) = 1 and t = ttk/n, for k > 1. 
If i = 1 and j = n then 

• mint f K -(l,n,t) — is attained when /ic n (l, 1,£) 
is maximum; 

• maxt /j*— (1, n, t) = 1 is attained when /ic n (l, l,t) 
is minimum. 

In all other cases, when i ^ j (and i ^ l,n) 

• mint f K -(i,j,t) = is attained when fK n (i,i,t) is 
maximum; 

• m&xt f K -(i, j,t) = 2/n is attained when fK n (i,i,t) 
is minimum. 

Proof. Let us define the n x n matrix P n such that 

[-Pn]l : l = [-Fn]n,n = lj [-Fn]l,n = = — lj an d 

[•Pn]ij = 0j otherwise. The Laplacian of K~ can be writ- 
ten as 

( n-2 -1 ••• -1 \ 
-1 ri-1-1 • • • -1 



HK-) = 



: -1 ■• -1 : 

-1 • • • -1 n- 1 -1 

\ -1 ••• -1 n-2/ 
= L{K n )-P n . 

The XYZ adjacency matrix of K~ has then the form 



14- 2i(Jf-) 



1 7„ - 2 (nJ» - J„ - P n ) 



H(K-) 



This matrix elements are of the form [H(K n )]i n = 
[H(K-)] n<1 = 0, [H(K-)] iti = m - 2d(i) ' and 
{H(K~)]ij — 2. The eigenvalues of H{K~) are then 

\f = n(n - l)/2 - 1, A™ = [H(K-)} 1A = n(n - 

l)/2 -271 + 3 and A,^ 2] = [H{K~)] 2 , 2 - [H{K-)\ lfl = 
n(n — 5)/2 — 1. We can chose an orthonormal basis of 
eigenvectors such that, from the spectral decomposition 



of the unitary matrix Ut(K n ) = e LH ( K n )*, we have the 
following diagonal entries: 



+ (| - l)e"M 
if i = 1, n: 



otherwise. 



Let us first consider the minimum value of f K - (i,i,t). 
In general, we would need to distinguish two cases de- 
pending on the parity of n. Here we take n to be 
a multiple of 4. This implies that n is always even. 
A simple calculation shows the next facts. It follows 
that mint f K - (1, 1, t) = and max t f K - (1, 1, t) = 1 if 
i = E + (Jfe > 0), and t = 7rfc/2 (fc > 1), respectively. 
The same holds for i — n. If i ^ l,n, it follows that 
min t /jj— (i, i, i) = 1 — 2/n and max t f K -(i,i,t) = 1 if 

i = j- + ^ and t = irk/n, respectively. 
The off-diagonal entries of U t {K~) are 



f n- 1 /V4©- 1 ]' - Hg-^KS)- 2 ^ 3 ]' 



[Ut{K-)]i 



+ (f -l) e -H— =- 
if i = 1 and j = n; 

-*[(;)-!]* _ 

, otherwise. 



Given this equation, for i = 1 and j = n, we 
have min t f K - (1, n, i) = when £ = 7rfc/2 and 

max t / Kr -(l,7i,t) = 1 when t = f + i.e., f K -(l,n,t) 
is minimum when /«— (1, 1, i) is maximum and mz. The 
last case to describe is when i 7^ j and i 7^ \,j 7^ n: 
min t j,t) = and max t f K - t) = 2/n if t = 

itk/n (k > 0) and t = ^ + ^ (fc > 1), respectively. ■ 

The fidelity given by the evolution in K~, with n = 
4,8,16,40 is illustrated in the figures below (top L/R, 
bottom L/R, resp.) for £ £ [0, tt]. The thick lines repre- 
sent f K - (1, 1, t); the thin ones, (1, n, i). The dashed 
lines represent f K - (2, 2, i); the dotted ones, f K - (2, 3, t). 
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FIG. 1: There is PST between vertex 1 and 8 in the top-left 
graph; for the pairs {1, 8} and {3, 6} in the top-right one; for 
the pairs {1, 8}, {2, 7}, {3, 6}, {4, 5} in the bottom graph. The 
dashed lines represent deleted edges. 



When n is a multiple of 4, the unitary Ut(K n ) can be 
used to route a qubit over n sites. To put in communica- 



tion site i and site j, we need to let the system evolve for 
a time t = \ + ^ , once we have deleted from K n exactly 
the edge if. 

The setting described above can be generalized by 
deleting more than a single edge from K n . In fact, as 
far as we delete edges without common vertices (or, in 
other words, vertex disjoint edges), we still can obtain 
PST between the end points of deleted edges, whenever 
the number of vertices in the graph is a multiple of 4. 
The maximum number of edges that can be deleted is 
then n/2. In this case, the deleted set corresponds to a 
perfect matching, that is, a set of vertex disjoint edges 
including all the vertices of the graph. As it is for K~, 
PST occurs when t = f + ^ ■ Figure Q] describes an 
example. 



V. CONCLUSION 

We have given an explicit relation between the Hamil- 
tonian for networks of interacting qubits in XYZ model 
and the graph Laplacian. We have studied evolution in 
networks modeled by the complete graph K n and the 
complete graph minus one edge K~ . While K n does not 
allow PST, we have seen that K ~ does allows PST be- 
tween the two non-adjacent sites, when n is a multiple of 
four. This result can be used to arbitrarily route a qubit 
between any two vertices of K~ . 
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